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This article presents the physics case for a new high-energy, ultra-high statistics neutrino scatter- 
ing experiment, NuSOnG (Neutrino Scattering on Glass). This experiment uses a Tevatron-based 
neutrino beam to obtain over an order of magnitude higher statistics than presently available for 
the purely weak processes + e~ — > z/ M + e~ and v^, + eT — + v a + . A sample of Deep In- 
elastic Scattering events which is over two orders of magnitude larger than past samples will also 
be obtained. As a result, NuSOnG will be unique among present and planned experiments for its 
ability to probe neutrino couplings to Beyond the Standard Model physics. Many Beyond Standard 
Model theories physics predict a rich hierarchy of TeV-scale new states that can correct neutrino 
cross-sections, through modifications of Zvv couplings, tree-level exchanges of new particles such 
as Z's, or through loop-level oblique corrections to gauge boson propagators. These corrections are 
generic in theories of extra dimensions, extended gauge symmetries, supersymmetry, and more. The 
sensitivity of NuSOnG to this new physics extends beyond 5 TeV mass scales. This article reviews 
these physics opportunities. 



I. INTRODUCTION 

Exploring for new physics at the "Terascale" - energy 
scales of ~ 1 TeV and beyond - is the highest priority 
for particle physics. A new, high energy, high statistics 
neutrino scattering experiment running at the Tevatron 
at Fermi National Accelerator Laboratory can look be- 
yond the Standard Model at Terascale energies by mak- 
ing precision electroweak measurements, direct searches 
for novel phenomena, and precision QCD studies. In this 
article we limit the discussion to precision electroweak 
measurements; QCD studies and their impact on the pre- 
cision measurements are explored in ref. [H [2] . The ideas 
developed in this article were proposed within the con- 
text of an expression of interest for a new neutrino ex- 



periment, NuSOnG (Neutrino Scattering On Glass) pQ. 

A unique and important measurement of the NuSOnG 
physics program is the ratio of neutral current (NC) and 
charged current (CC) neutrino-electron scattering, which 
probes new physics. The leading order Feynman dia- 
grams for these processes are shown in Fig. [l] The NC 
process, + er — > + e~ , called "elastic scattering" 
or ES, provides the sensitivity to the Terascale physics. 
This process can explore new physics signatures in the 
neutrino sector which are not open to other, presently 
planned experiments. The CC process, called "inverse 
muon decay" or IMD, + e~ — ► v e + pT , is well un- 
derstood in the Standard Model due to precision mea- 
surement of muon decay [3] . Since the data samples are 
collected with the same beam, target and detector at 



FIG. 1: Left: "elastic scattering" (ES). Right: "Inverse Muon 
Decay" (IMD). 
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the same time, the ratio of ES to IMD events cancels 
many systematic errors while maintaining a strong sen- 
sitivity to the physics of interest. Our measurement goal 
of the ES to IMD ratio is a 0.7% error, adding system- 
atic and statistical errors in quadrature. The high sen- 
sitivity which we propose arises from the combined high 
energy and high intensity of the NuSOnG design, leading 
to event samples more than an order of magnitude higher 
than past experiments. 

Normalizing the ES to the IMD events represents an 
important step forward from past ES measurements, 
which have normalized neutrino-mode ES measurements 
to the antineutrino mode, P M + e~ — > + e~[4, 5]. The 
improvement is in both the experimental and the theo- 
retical aspects of the measurement. First, the flux con- 
tributing to IMD and v ES is identical, whereas neutrino 
and antineutrino fluxes are never identical and so require 
corrections. Second, the ratio of v ES to v ES cancels 
sensitivity to Beyond Standard Model (BSM) physics ef- 
fects from the NC to CC coupling ratio, p, which are 
among the primary physics goals of the NuSOnG mea- 
surement. In contrast, there is no such cancellation in 
the ES to IMD ratio. 

The design of this experiment, described in Sec. [IT] 
is driven both by requiring sufficient statistics to make 
precision neutrino-electron scattering measurements and 
by the need for a neutrino flux which does not extend 
below the IMD threshold. The threshold for IMD events 
is 
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E„>E ll >-^ = 10.9 GeV, (1) 
Zm e 

where we have dropped the small m 2 term for simplicity. 
The functional form above threshold, shown in Fig. [2j is 
given by (1 — m 2 ^/ E 2 m ) 2 , where E cm is the center of mass 
energy. Thus a high energy neutrino beam is required to 
obtain a high statistics sample of these events. The flux 
design should provide a lower limit on the beam energy 
of about 30 GeV, still well above the IMD threshold. 

Sec. |III| describes the Standard Model Physics of neu- 
trino electroweak scattering, for both electron and quark 
targets. In this section, the value of the normalization of 
the ES to IMD events is further explored. The very high 
statistics will also permit an electroweak measurement 
using the deep inelastic scattering (DIS) data sample 
from NuSOnG, via the "Paschos Wolfenstein method" 



FIG. 2: Threshold factor for the IMD cross section, as a 
function of neutrino energy. 

(PW) [Qj. The best electroweak measurement using 
DIS events to date comes from the NuTeV experiment, 
which has observed an anomaly. The status of this re- 
sult is reviewed below. Making conservative assumptions 
concerning systematic improvements over NuTeV, our 
measurement goal using this technique is a 0.4% error 
on sin 2 #vKi adding statistical and systematic errors in 
quadrature. 

In Sec.|lVj we discuss NuSOnG 's potential to discover 
or constrain new physics through indirect probes, by 
making precision measurements of SM processes to look 
for deviations from SM predictions. We first frame the is- 
sue by considering in turn several model-independent pa- 
rameterizations of possible new physics and asking what 
constraints will be imposed on new physics in the event 
NuSOnG agrees with the SM. (1) Oblique correction 
parameters describe the effects of heavy new states in 
vector boson loops. (2) New states may induce higher- 
dimensional effective operators involving neutrinos. Fi- 
nally, (3) new states may modify the couplings of the 
gauge bosons to neutrinos and leptons, including possi- 
bly violating lepton universality. In each case we consider 
the ability of NuSOnG to detect or constrain these types 
of deviations from the SM. 

In Sec. |Vj we examine specific models for new physics. 
We begin by presenting the sensitivity to a set of new 
physics models. In particular, we consider 

• typical Z' models, 

• non-degenerate leptoquark models, 

• R-parity violating SUSY models, 

• extended Higgs models. 

The models were selected because they are often used as 
benchmarks in the literature. While this list is not ex- 
haustive, it serves to illustrate the possibilities. For each 
case, we consider how NuSOnG compares to other mea- 
surements and note the unique contributions. We end 
this section by approaching the question from the oppo- 
site view, asking: how could the results from NuSOnG 
clarify the underlying physics model, should evidence of 
new physics emerge from LHC in the near future? 
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FIG. 3: The assumed energy-weighted flux, from the NuTeV 
Experiment [7 , in neutrino mode (left) and antineutrino 
mode (right). Black: muon neutrino, red: muon antineutrino, 
blue: electron neutrino and antineutrino flux. 



Two further studies which can be performed by Nu- 
SOnG are QCD measurements and direct searches. The 
very large (~ 600 million event) DIS sample will allow 
the opportunity for precision studies of QCD. There are 
many interesting measurements which can be made in 
their own right and which are important to NuSOnG 's 
Terascale physics program. The very high flux will also 
permits direct searches for new physics. Those which 
complement the physics discussed in this paper include: 

• non-unitarity in the light neutrino mixing matrix; 

• wrong-sign inverse muon decay (WSIMD), + 
e~ -> M~ + 

• decays of neutrissimos, i.e., moderately-heavy 
neutral-heavy-leptons, with masses above 45 GeV. 

For more information on these studies, see refs. [HH]. 



II. CONCEPTUAL DESIGN FOR THE 
EXPERIMENT 

In order to discuss the physics case for a new high en- 
ergy, high statistics experiment, one must specify certain 
design parameters for the beam and detector. The beam 
and detector should marry the best aspects of NuTeV 
[7], the highest energy neutrino experiment, and Charm 
II [3], the experiment with the largest ES sample to date. 
The plan presented here is not optimized, but provides a 
basis for discussion. The final design of the NuSOnG de- 
tector will be based on these concepts, and is still under 
development . 

In this section, we present, but do not justify, the de- 
sign choices. Later in this article, we discuss the reason- 
ing for the choices, particularly in Sees . |III C] and |TlI D | 

We will assume a beam design based on the one used by 
the NuTeV experiment [7] , which is the most recent high 
energy neutrino experiment. This experiment used 800 
GeV protons on target. The beam flux, shown in Fig. [3] 
is ideal for the physics case for several reasons. There 
is essentially no flux below 30 GeV, hence all neutrinos 



are well above the IMD threshold. It is sign-selected: in 
neutrino mode, 98.2% of neutrino interactions were due 
to n + and K + secondaries, while in antineutrino mode 
97.3% came from ir~ and K~. The "wrong sign" content 
was very low, with a 0.03% antineutrino contamination 
in neutrino mode and 0.4% neutrino contamination in 
antineutrino mode. The electron-flavor content was 1.8% 
in neutrino mode and 2.3% in antineutrino mode. The 
major source of these neutrinos is K e3 decay, representing 
1.7% of the total flux in neutrino mode, and 1.6% in 
antineutrino mode. 

Redesign of the beamline for NuSOnG is expected to 
lead to modest changes in these ratios. For example, 
if the decay pipe length is 1.5 km rather than 440 m, 
as in NuTeV, the tt/K ratio increases by 20% and the 
fractional v e content is reduced. 

With respect to Tevatron running conditions, we will 
assume that twenty times more protons on target (POT) 
per year can be produced for NuSOnG compared to 
NuTeV. This is achieved through three times higher in- 
tensity per pulse (or "ping"). Nearly an order of mag- 
nitude more pulses per spill are provided. Our studies 
assume 4 x 10 19 POT/year, with 5 years of running. 
Preliminary studies supporting these goals are provided 
in ref. [8]. 

The event rates quoted below are consistent with 
1.5xl0 20 protons on target in neutrino running and 
0.5 x 10 20 protons on target in antineutrino running. The 
choice to emphasize neutrino running is driven by obtain- 
ing high statistics ES, which has a higher cross section for 
neutrino scatters, and to use the IMD for normalization - 
this process only occurs in neutrino scattering. The Stan- 
dard Model forbids an IMD signal in antineutrino mode. 
However, some antineutrino running is required for the 
physics described in the following sections, especially the 
PW electroweak measurement. 

The beam from such a design is highly forward di- 
rected. NuTeV was designed so that 90% of the neutrinos 
from pion decay were contained within the detector face, 
where the detector was located at 1 km. For NuSOnG, 
which will use a 5 m detector, ~90% of the neutrinos 
from pion decay are contained at ^3 km. 

The optimal detector is a fine-grained calorimeter 
for electromagnetic shower reconstruction followed by 
a toroid muon spectrometer. This allows excellent re- 
construction of the energy of the outgoing lepton from 
charged current events. We employ a Charm II style de- 
sign [S], which uses a glass target calorimeter followed 
by a toroid. We assume one inch glass panels with ac- 
tive detectors interspersed for energy and position mea- 
surement. Glass provides an optimal choice of density, 
low enough to allow electromagnetic showers to be well 
sampled, but high enough that the detector length does 
not compromise acceptance for large angle muons by the 
toroid. Approximately 10% of the glass will be doped 
with scintillator to allow for background studies, as dis- 
cussed in Sec. IIII Dl 

The design introduces four identical sub-detectors of 
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Vfj, CC Deep Inelastic Scattering 
NC Deep Inelastic Scattering 
v^, electron NC elastic scatters (ES) 
Vp electron CC quasi-elastic scatters (IMD) 
Vfx CC Deep Inelastic Scattering 
Tip, NC Deep Inelastic Scattering 
electron NC elastic scatters (ES) 
Up electron CC quasi-elastic scatters (WSIMD) 



and fermion / leads to the effective interaction: 



TABLE I: Rates assumed for this paper. NC indicates "neu- 
tral current" and CC indicates "charged current." 



this glass-calorimeter and toroid design, each a total of 
29 m in length (including the toroid) . Between each sub- 
detector is a 15 m decay region for direct searches for 
new physics. The total fiducial volume is 3 ktons. 

The NuSOnG run plan, for reasons discussed in 
Sec " 



IIIB and IIIC 



run plan, tor reasons 
concentrates on running in neutrino 



mode. This design will yield the rates shown in Table [T] 
These rates, before cuts, are assumed throughout the rest 
of the discussion. We can compare this sample to past 
experiments. The present highest statistics sample for 
and ES is from CHARM II, with 2677±82 events in 
neutrino mode and 2752±88 events in antincutrino mode 
[5]. Thus the proposed experiment will have a factor of 
30 (2.5) more i/(j>)-electron events. As an example, after 
cuts, the first method of analysis described in Sec. 
retains 63% of the v sample. 



For deep inelastic scatter- 
ing, 600M and 190M events are expected in neutrino and 
antincutrino modes, respectively. After minimal cuts to 
isolate DIS events [TO], NuTeV had 1.62M DIS (NC+CC) 
events in neutrino mode and 0.35M in antineutrino mode; 
thus NuSOnG has orders of magnitude more events. 

The detector will incorporate several specialized re- 
gions. A region of fine vertex-tracking facilitates mea- 
surements of the strange sea relevant for the electroweak 
analysis, as described in ref. [2J. Two possibilities are 
under consideration: an emulsion detector or a silicon 
detector of the style of NOMAD-STAR [11 . Both are 
compact and easily accommodated. For further QCD 
studies, it will also be useful to intersperse alternative 
target materials: C, Al, Fe, and Pb [2J. 



III. ELECTROWEAK MEASUREMENTS IN 
NEUTRINO SCATTERING 

Neutrino neutral current (NC) scattering is an ideal 
probe for new physics. An experiment like NuSOnG is 
unique in its ability to test the NC couplings by studying 
scattering of neutrinos from both electrons and quarks. A 
deviation from the Standard Model predictions in both 
the electron and quark measurements would present a 
compelling case for new physics. 

The exchange of the Z boson between the neutrino v 
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where the Standard Model values of the couplings are: 
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or equivalcntly, 
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Here, l( and Qf are the weak isospin and electromagnetic 
charge of fermion /, respectively. In these formulas, p is 
the relative coupling strength of the neutral to charged 
current interactions (p = 1 at tree level in the Standard 
Model). The weak mixing parameter, sin 2 9yy, is related 
(at tree level) to Gp , Mz and a by 

47TCV 

V2G F M 2 ' 



sin 2 26 w = 



(5) 



A. Neutrino Electron Elastic Scattering 

The differential cross section for v„ and P M ES, defined 
using the coupling constants described above, is: 



do 
dT 



2G%m e 
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+{9l9vT9l9 A f 
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m e T 



(6) 



The upper and lower signs correspond to the neutrino 
and anti-neutrino cases, respectively. In this equation, 
E v is the incident energy and T is the electron recoil 
kinetic energy. 

More often in the literature, the cross section is defined 
in terms of the parameters (gy e ,gA e )i which are defined 
as 
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In terms of these parameters, we can write: 
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When m e <C E u , as is the case in NuSOnG, the third 
term in these expressions can be neglected. If we intro- 
duce the variable y = T / E Vl then 



da G 2 F m e E v 
dy 2ir 



(g^±g^f + (g^T9 V A e f^-v) 2 



(9) 



Integrating, we obtain the total cross sections which are 
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Note that 
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The ratio of the integrated cross sections for neutrino 
to antineutrino electron ES is 



a{y^ e) 
Rv/v = — ~ — r = 3 



l-4sin 2 W + ^sin 4 9 W 



l-4sin 2 
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(13) 



Fig. [ijtop) shows the results for sin 2 9w from many past 
experiments which have used this a v/v ES ratio." 

In the ratio, R v /a, the dependence on p canceled. This 
directly extracts sm 2 9w- The relationship between the 
error on the ratio and the error on sin 2 0^/, which for 
convenience we abbreviate as z, is: 
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FIG. 4: Measurements of sin 2 8w from past experiments. 
Top: neutrino-electron elastic scattering experiments. Bot- 
tom: neutrino DIS experiments. All DIS results are adjusted 
to the same charm mass (relevant for experiments not using 
the PW method). The Standard Model value, indicated by 
the line, is 0.2227 [12]. 



for z = 0.2227 (or R u/P = 1.242). Roughly, the fractional 
error on sin 2 9w is 60% of the fractional error on R„/p. 



B. A New Technique: Normalization Through 
IMD 



An experiment such as NuSOnG can make indepen- 
dent measurements of the electroweak parameters for 
both V/j, and i/^-electron scattering. We can achieve this 
via ratios or by direct extraction of the cross section. In 
the case of ^-electron scattering, we will use the ratio of 
the number of events in neutrino-electron elastic scatter- 
ing to inverse muon decay: 



U NC 
T IMD 



X 
X 



(16) 



Because the cross section for IMD events is well deter- 
mined by the Standard Model, this ratio should have low 
errors and will isolate the EW parameters from NC scat- 
tering. In the discussion below, we will assume that the 
systematic error on this ratio is 0.5%. 

In the case of data, the absolute normalization 
is more complex because there is no equivalent process 
to inverse muon decay (since there are no positrons in 
the detector). One can use the fact that, for low ex- 
change energy (or "nu") in Deep Inelastic Scattering, 
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FIG. 5: Kinematic distributions for IMD events from incident neutrino energy between 100 and 200 GeV. Left: y distribution; 
right: 9u distribution. Black: distribution of events before cuts; Red: distribution after cuts for analysis method 1 (see 
SeclmDi. 



the cross sections in neutrino and antincutrino scatter- 
ing approach the same constant, A [13] . This is called 
the "low nu method" of flux extractions. For DIS events 
with low energy transfer and hence low hadronic en- 

: <b u A and 



ergy (5 < E had < 10 GeV), 

N l °™ I g had = <S>»A. The result is that the electroweak 
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parameters can be extracted using the ratio 



j.rlow E h 

^uDIS 



-) 



*jlow E ha 
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V v e ) 



JVC 



X $ p 



a IMD x 

(17) 

The first ratio cancels the DIS cross section, leaving the 
energy-integrated v to v flux ratio. The IMD events in 
the denominator of the second term cancel the integrated 
v flux. The NC elastic events cancel the integrated v flux. 

Because of the added layer of complexity, the antineu- 
trino ES measurement would have a higher systematic 
error than the neutrino ES scattering measurement. The 
potentially higher error is one factor leading to the plan 
that NuSOnG concentrate on neutrino running for the 
ES studies. 

As shown in Fig. [2j IMD events have a kinematic 
threshold at 10.9 GeV. These events also have other in- 
teresting kinematic properties. The minimum energy of 
the outgoing muon in the lab frame is given by 



Timlin 
lab 



2m P 



10.9 GeV. 



(18) 



In the detector described above, muons of this energy 
and higher will reach the toroid spectrometer without 
ranging-out in the glass. An interesting consequence is 
that, independent of E v , the energy transfer in the inter- 
action has a maximum value of 



Vr, 



1 



10.9 GeV 



(19) 



Thus at low the cutoff in y is less than unity, as 
shown in Fig. [5] (left). The direct consequence of this is 
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FIG. 6: Reconstructed neutrino energy (red) for IMD events 
before cuts compared to true neutrino energy (black). 



a strong cutoff in angle of the outgoing muon, shown in 
Fig. [5] (right). In principle, one can reconstruct the full 
neutrino energy in these events: 



E 



IMD 



1 2m e E a -ml - ml 



2 m. - E„ 



p M cos 9^ 



(20) 



This formula depends on 9^, which is small. The recon- 
structed E v is smeared by resolution effects as seen in 
Fig. [6] While the analysis can be done by summing over 
all energies, these distributions indicate that an energy 
binned analysis may be possible. This is more powerful 
because one can fit for the energy dependence of back- 
grounds. For the illustrative analyzes below, however, we 
do not employ this technique. 

The error on sin 2 6w extracted from this ratio, 
Res/imdi assuming a Standard Model value for p, is 
the same as the error on the ratio: 



5(sin 2 9 w ) ^ 5R ES / IMD 
sin 2 9 w Res/imd 



(21) 
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Ref. [2] provides a useful summary of radiative cor- 
rections for the ES and IMD processes, which were orig- 
inally calculated in Ref. [15]. The error from radiative 
corrections is expected to be below 0.1%. It is noted 
that to reduce the error below 0.1%, leading two-loop ef- 
fects must be included. A new evaluation of the radiative 
corrections is underway |16j . 



C. IMD Normalization vs. v Normalization 

NuSOnG can measure both the vjv ES ratio, as in the 
case of past experiments shown in Eq. ( 13 ), as well as the 



ES/IMD ratio. In the case of the former, to obtain the 
best measurement in a 5 year run, one would choose a 1:3 
ratio of run time in v versus v mode. In the latter case, 
one would maximize running in v mode. The result of the 
two cases is a nearly equal error on sin 2 6\y , despite the 
fact that the error on the vjv ES is nearly twice that of 
the ES/IMD ratio. To understand this, compare Eq. ( 15 ) 
to Eq. ( |2l| ). However, the ES/IMD ratio is substantially 
stronger for reasons of physics. Therefore, our conceptual 
design calls for running mainly with a v beam. In this 
section we explore the issues for these two methods of 
measurement further. We also justify why the precision 
measurement requires high energies, only available from 
a Tevatron-based beam. 



1. Comparison of the Two Measurement Options 

From the point of view of physics, The ES/IMD ra- 
tio is more interesting than the vjv ES ratio. This is 
because p has canceled in the vjv ES ratio of Eq. 
leaving the ratio insensitive to physics which manifests 
itself through changes in the NC coupling. Many of the 
unique physics goals of NuSOnG, discussed in Sec. |IV| 
depend upon sensitivity to the NC coupling. 

An equally important concern was one of systemat- 
ics. The v and v fluxes for a conventional neutrino beam 
are substantially different. For the case of NuSOnG, the 
fluxes are compared in Fig. [3] Predicting the differences 
in these fluxes from secondary production measurements 
and simulations leads to substantial systematic errors. 
For beams at high energies (> 30 GeV), such as Nu- 
SOnG, the "low nu" method [13] for determining the 
ratio of the neutrino to antineutrino fluxes from Deep 
Inelastic events, developed by CCFR and NuTeV and 
described in Sec. |IIIB| can be employed. However, this 
leads to the criticism that one has introduced a new pro- 
cess into the purely-leptonic analysis. 

Neither criticism is relevant to the ES/IMD ratio. The 
sensitivity to the new physics through the couplings does 
not cancel. Because both processes are in neutrino mode, 
the flux exactly cancels, as long as the neutrino energies 
are well above the IMD thresho ld (thi s will be illustrated 
in the analysis presented in Sec. HID ). This ratio has the 



added advantage of needing only neutrino- mode running, 



which means that very high statistics can be obtained. 
This is clearly the more elegant solution. 

It should be noted that nothing precludes continued 
running of NuSOnG beyond the 5-year plan presented 
here. This run-length was selected as "reasonable" for 
first results. If interesting physics is observed in this first 
phase, an extended run in antineutrino mode may be 
warranted, in which case both the ES/IMD and vjv ES 
ratios could be measured. The latter would then con- 
strain sin 2 9w in a pure neutrino measurement and the 
former is then used to extract p. 

To measure the ES/IMD ratio to high precision, there 
must be little low energy flux. This is because the IMD 
has a threshold of 10.9 GeV, and does not have substan- 
tial rate until ~ 30 GeV. The low-energy cut-off in the 
flux (see Fig. [3) coming from the energy-angle correlation 
of neutrinos from pion decay, is ideal. 



2. Why a Tevatron-based Beam is Best for Both Options 

The ES/IMD measurement is not an option for the 
planned beams from the Main Injector at Fermilab. For 
both presently planned Main Injector experiments at Fer- 
milab T7] and for the proposed Project-X DUSEL beam 
[IS] , the neutrino flux is peaked at ^5 GeV. The major- 
ity of the flux of these beams is below 5 GeV, and most 
of the flux is below the 10.9 GeV IMD threshold. Be- 
cause of this, one simply cannot use the IMD events to 
normalize. 

In principle, the vjv ES ratio could be used. How- 
ever, in practice this will have large systematics. The v 
and v fluxes for a horn beam are significantly different. 
First principles predictions of secondary mesons are not 
sufficient to reduce this error to the precision level. The 
energy range is well below the deep inelastic region where 
the "low nu" method can be applied to accurately extract 
a vjv flux ratio. Other processes, such as charged-current 
quasi-elastic scattering, could be considered for normal- 
ization, but the differences in nuclear effects in neutrino 
and antineutrino scattering for these events is not suffi- 
ciently well understood to yield a precision measurement. 

Lastly, the ES rates for the present Main Injector 
beams are too low for a high statistics measurement. 
This is because the cross section falls linearly with en- 
ergy. Event samples on the order of 10k may be possible 
with extended running in the Project X DUSEL beam 
in the future. From the point of view of statistics, even 
though two orders of magnitude more protons on target 
are supplied in such a beam, the Tevatron provides a 
substantially higher rate of ES per year of running. 

Compared to the Main Injector beam, a Tevatron- 
based beam does not face these issues. The choice of 
running in neutrino mode provides the highest precision 
measurement while optimizing the physics. 
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Quantity 


Assumed Value 


Uncertainty 


Source of Estimate 


Muon 








Energy Resolution 


SE/E = 10% 


2.5% 


NuTeV testbeam measurement 


Energy Scale Error 




0.5% 


NuTeV testbeam measurement 


Angular Resolution 


59 = 0.011/S 96 rad 


2.5% 


Multiple scattering fit simulation 


Electron 








Energy Resolution 


5E/E = 0.23/E 05 


1.0% 


Same as CHARM II 


Energy Scale Error 


vec 1.0 X -E~t TUG 


1.0% 


Scaled from CHARM II with NuSOnG statistics 


Angular Resolution 


86 = 0.008/S ' 5 rad 


2.5% 


2 better than CHARM II due to sampling 


Flux 








Normalization 


1.0 


3% 


Current total cross section uncertainty 


Shape Uncertainty 


1.0 


1% 


Similar to NuTeV low-nu method 


Backgrounds 








CCQE 


1.0 


5% 


Extrapolated from NuTeV 


v s CCQE 


1.0 


3% 


Extrapolated from CHARM II 



TABLE II: Resolutions and systematic uncertainty estimates used in the parameterized Monte Carlo studies. The NuTeV 
estimates are based on Ref. |19| and the CHARM II estimates from Ref. Units for angles are radians and energies are in 
GeV. 



D. A 0.7% Measurement Goal for the ES to IMD 
Ratio 

Achieving 0.7% precision on the ES/IMD measure- 
ment depends on reducing the backgrounds to an accept- 
able level without introducing significant systematics and 
while maintaining high signal statistics. Many of the 
systematic uncertainties will tend to cancel. The most 
important background for both the v-e neutral current 
and IMD events comes from charged current quasi-elastic 
(CCQE) scatters {y e n — > pe and v^n — * pfi). These 
background CCQE processes have a much broader Q 2 as 
compared to the signal processes and, therefore, can be 
partially eliminated by kinematic cuts on the outgoing 
muon or electron. Initial cuts on the scattering angle 
and energy of the outgoing muon or electron can easily 
reduce the CCQE background by factors of 60 and 14 
respectively while retaining over 50% of the v-e neutral 
current and IMD signal. This leaves events with very 
forward scatters and outgoing scattered protons of low 
kinetic energy. 

Because the NuSOnG design is at the conceptual stage 
and in order to be conservative, we have developed two 
different strategies for achieving a 0.7% error. This serves 
as a proof of principle that this level of error, or better, 
can be reached. The first method relies on detecting pro- 
tons from the quasi-elastic scatter. The second method 
uses the beam kinematics to cut the low energy flux which 
reduces the CCQE background. 

These methods were checked via two, independently 
written, parameterized Monte Carlos. The parameter- 
ized Monte Carlos made the assumptions given in Table 
HTl where both the assumed values and uncertainties are 
presented. These estimates of resolutions and system- 
atic errors are based on previous experimental measure- 
ments or on fits to simulated data. One Monte Carlo 
used the Nuance event generator [20 to produce events, 
while the other was an independently written event gen- 



erator. Both Monte Carlos include nuclear absorption 
and binding effects. 

The first strategy uses the number of protons which 
exit the glass to constrain the total rate of the back- 
ground. In ~ 33% of the events, a proton will exit the 
glass, enter a chamber and traverse the gas. This sam- 
ples protons of all energies and Q 2 , since the interactions 
occur uniformly throughout the glass. After initial cuts, 
the protons are below 100 MeV, and therefore highly ion- 
izing. If we define 1 MIP as the energy deposited by a 
single minimum ionizing particle, like a muon, then the 
protons consistently deposit greater than 5 MIPs in the 
chamber. Thus, one can identify CCQE events by re- 
quiring >4 MIPS in the first chamber. The amount of 
remaining CCQE background after this requirement can 
be measured if a fraction such as 10% of the detector is 
made from scintillating glass that can directly identify 
CCQE events from light associated with the outgoing 
proton. A wide range of scintillating glasses have been 
developed [5T] for nuclear experiments. These glasses are 
not commonly used in high energy physics experiments 
because the scintillation time constant is typically on the 
order of 100 ns. In a neutrino experiment, which has 
inherently lower rates than most particle experiments, 
this is not an issue. CCQE events can be identified by 
the scintillation light from the proton assuming reason- 
able parameters for the glass and readout photomulti- 
plier tubes: 450 photons/MeV, an attenuation length of 
2 m, eight phototubes per glass sheet, quantum efficiency 
of the tubes of 20%. Using the identified CCQE events 
from the instrumented glass, the uncertainty in the resid- 
ual background can be reduced to 2.0% for the IMD mea- 
surement. For the CCQE background to the v^-e neutral 
current measurement, the uncertainty is assumed to be 
3% for the Monte Carlo prediction. Combining all the 
systematic errors leads to a ^0.7% accuracy on the v-e 
measurement as shown in Tab. IIIII 



In Tab. Ill the cancellation of the flux errors should be 
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noted. This occurred because we use the ES/IMD ratio, 
as discussed in the previous section. 

The second strategy involves reducing the relative 
CCQE background to signal by using a harder flux for the 
analysis. This study used the same Monte Carlos, with 
the resolutions listed in Tab. |TTJ as the first analysis. The 
total systematic and statistical error achieved was 0.6%. 
Below, we explain how a harder flux is obtained for the 
analysis. Then, we explain how this flux improves the 
signal-to-background in both the ES and IMD analyzes. 

The strong correlation between energy and angle at 
the NuSOnG detector is used to isolate the harder flux. 
This is simplest to express in the non-bend view of the 
beamline, where it is given for pions by the well-known 
off-axis formula: 



E„ 



0.43K 
1 + 7 2 (9 



(22) 



where 9 is the off-axis angle, 7 = E^/m^, is the en- 
ergy of the pion and E v is the energy of the neutrino. For 
the NuTeV beam and detector lay-out, this angle-energy 
dependence resulted in the sharp cutoff of the flux for 
< 30 GeV shown in Fig. [3j Using the NuTeV G3 beam 
Monte Carlo [7], we have shown that by selecting ver- 
tices in the central region of the detector, one can adjust 
the energy where the flux sharply cuts off. Adjusting the 
aperture to retain flux above 50 GeV reduces the total 
event rate by 55%. 

A harder flux allows for background reduction in both 
the ES and the IMD samples while maintaining the sig- 
nal at high efficiency. In the case of ES events, the back- 
ground is from v e CCQE. The energy distribution of the 
electron is substantially different in the two cases. In the 
case of v e CCQE events, the electron carries most of the 
energy of the incoming neutrino because the exchange en- 
ergy in the interaction is small. Thus the CCQE events 
produced by the harder flux populate the visible energy 
range above 50 GeV. On the other hand, the outgoing 
electron in ES events tends to populate the low visible 
energy region due to the combination of a flat y distribu- 
tion for the process convoluted with the incident neutrino 
energy spectrum. The result is that a cut on the visible 
energy less than 50 GeV reduces the error from the v e 
CCQE background to a negligible level. To understand 
the improvement in the IMD analysis, consider Fig. [2] 
which shows the threshold effects. The IMD signal is also 
rising with energy In contrast, the CCQE rate, which 
is the most significant background, is flat with energy for 
fluxes above 1 GeV. This signal-to-background is greatly 
improved with a high energy flux. This allows looser cuts 
to be applied, which in turn reduces the systematics. 

These two analyzes use substantially different strate- 
gies and can, in principle, be combined. Given these 
preliminary studies, we feel confident that as the detec- 
tor moves from a conceptual to real design, we will be 
able to achieve a better than 0.7% error. However, for 
this paper we take the conservative approach of assuming 
0.7%. 



E. Neutrino Quark Scattering 



Substantially higher precision has been obtained us- 
ing neutrino-quark scattering, which compares neutral- 
current (NC) to charged-current (CC) scattering to ex- 
tract sin 2 9\y . However, these experiments are subject to 
issues of modeling in the quark sector. Fig. ^bottom) 
reviews the history of these measurements. 

The lowest systematic errors come from implement- 
ing a "Paschos-Wolfenstein style" [B] analysis. This PW 
technique would be used by any future experiment, in- 
cluding NuSOnG. This requires high purity v and v 
beams, for which the following ratios of DIS events could 
be formed: 



R u 



R? 



'NC 



'CC 



'NC 
T CC 



Paschos and Wolfenstein [6| recast these as: 



R- = 



'NC 



'NC 



CC 



R v - rR* 
1 -r : 



(23) 
(24) 

(25) 



where r = G v cc j a v cc . In Rr many systematics cancel 
to first order, including the effects of the quark and an- 
tiquark seas for u, d, s, and c. Charm production only 
enters through d va i ence (which is Cabibbo suppressed) 
and at high x; thus the error from the charm mass is 
greatly reduced. The cross section ratios can be written 
in terms of the effective neutrino-quark coupling param- 
eters g\ and g\ as 



in which 



R" 
It* 

R- 



A 



A 



A 
A 



A 



r 9R 
1 2 

~9r 



■At 



P% 



(26) 
(27) 

(28) 



(2<7l5l) 2 
{2glglf + {2glg d R f 

P (g Sln VW)- 



5 . 4 

sin i 



■6 



w 



W I 



(29) 
(30) 



In a variation on the PW idea, rather than directly 
form R~, NuTeV fit simultaneously for R v and R v to 
extract sin 2 6w, obtaining the value sin 2 9yy = 0.2277 ± 
0.00162. Events were classified according to the length 
of hits in the scintillator planes of the NuTeV detector, 
with long events identified as CC interactions and short 
events as NC. An important background in the CC sam- 
ple came from pion decay-in-flight, producing a muon in 
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TMD TInpprtaintv 

J..1V1-1 ' KJ LLV^VjL L/CLlllU y 


T-^S TInpprtaintv 


TInpprtaintv on Ratio 


Statistical Uncertainty 


0.18% 


0.46% 


0.49% 


R pcnl n + i nn ^tti pa n n cr 

l_\,COLJl LI L-1U11 kJllACCll lilt. 








^Cp i _ +o 5% 

(_/ ^ J_J^ y — . U /(J 


0.00% 


0.00% 


0.00% 


S(0 M ) = ±2.5% 


0.04% 


0.00% 


0.04% 




0.00% 


0.01% 


0.01% 


1 h2 5% 


0.00% 


0.09% 


0.09% 


Knormr s,/ 1 ^ lo 
Ljllclgy OLdlL 








<5(Escale M ) = 0.5% 


0.37% 


0.00% 


0.37% 


rSfEsoale ) — 1 5% 


0.00% 


0.19% 


0.19% 


Flux 








ISTnrTTifil i 7 a tinn 

1 > VJ1111c1j11ZjCLLi1W11 


3.00% 


3.00% 


0.00% 


High energy flux up 1% 


0.25% 


0.25% 


0.00% 


Low energy flux up 1% 


0.15% 


0.13% 


0.02% 


IMD Background: statistical error 


0.06% 


0.00% 


0.06% 


2.0% systematic error 


0.26% 


0.00% 


0.26% 


f M e Background: statistical error 


0.00% 


0.12% 


0.12% 


3% systematic error 


0.00% 


0.19% 


0.19% 




Total Syst. 


Uncertainty on Ratio 


0.54% 




Total Stat. 


Uncertainty on Ratio 


0.51% 




Total Uncertainty on Ratio 


0.74% 



TABLE III: Estimates of the IMD and ES uncertainties using a > 5 MIP cut on the first downstream chamber. The columns 
give the errors for each process and then for the ratio. Errors are included for statistical uncertainties and uncertainties 
associated with the knowledge of resolution smearing, energy scale, flux shape, and backgrounds. The flux shape uncertainties 
are significantly reduced in the ratio measurement. 



a NC shower. Significant backgrounds in the NC sample 
came from muons which ranged out or exited and from 
u e CC scatters which do not have a muon and thus are 
classified as "short." 

In this paper, we present the sensitivity of NuSOnG 
to new physics if the NuTeV errors are reduced by a 
factor of ~ 2. This is a very conservative estimate, since 
most of the improvement comes from higher statistics. 
Only a 90% improvement in the systematics is required 



to reach this goal. Tab. IV argues why a 90% reduction 
in systematic error should be straightfroward to achieve. 
It is likely that the NuSOnG errors will be lower, but this 
requires detailed study. 

In Table |IV[ we list the errors which NuTeV identi- 
fied in their original analysis and indicate how NuSOnG 
will improve each error. Many of the largest experimen- 
tal systematics of NuTeV are improved by introducing a 
fine-grained sampling calorimeter. The NuTeV detector 
had four inches of iron between unsegmented scintillator 
planes and eight inches between drift chamber planes. 
Better lateral segmentation and transverse detection will 
improve identification of scatters from intrinsic v e s in the 
beam and separation of CC and NC events by improved 
three-dimensional shower shape analyzes. The NuTeV 
analyzes of the intrinsic v e content [22] and the CC/NC 
separation for the sin 2 Q\y analysis which relied strictly 
on event length. With this said, the power of classifying 
by event length is shown by the fact that the NuTeV in- 
trinsic v e analysis was sensitive to a discrepancy in the 
predicted intrinsic v e rate which was recently resolved 
with a new measurement of the K e % branching ratio that 



was published in 2003. Details of these issues are consid- 
ered in the next section. 



F. The NuTeV Anomaly 

From Fig. [4j it is apparent that the NuTeV measure- 
ment is in agreement with past neutrino scattering re- 
sults, although these have much larger errors; however, in 
disagreement with the global fits to the electroweak data 
which give a Standard Model value of sin 2 6\y = 0.2227 
[251. Expressed in terms of the couplings, NuTeV mea- 
sures: 



0.30005 ±0.00137 



g 2 R = 0.03076 ±0.00110, 



(31) 
(32) 



which can be compared to the Standard Model values of 
g\ = 0.3042 and g\ = 0.0301, respectively. 

NuTeV is one of a set of Q 2 -C m? z experiments mea- 
suring sin 2 6w- It was performed at Q 2 — 1 to 140 GeV 2 , 
(Ql) = 26 GeV 2 , (Q 2 ) = 15 GeV 2 , which is also the 
expected range for NuSOnG. Two other precision low 
Q 2 measurements are from atomic parity violation |26j 
(APV), which samples Q 2 - 0; and SLAC E158, a M0ller 
scattering experiment at average Q 2 — 0.026 GeV 2 |27j . 
Using the measurements at the Z-pole with Q 2 = M 2 
to fix the value of sin 2 ^, and evolving to low 2 |28|. 
the APV and SLAC E158 are in agreement with the 
Standard Model. However, the radiative corrections to 
neutrino interactions allow sensitivity to high-mass par- 
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Source 


NuTeV 
Error 


Method of reduction in NuSOnG 


Statistics 


0.00135 


Higher statistics 


Ve, flux prediction 


0.00039 


Improves in-situ measurement of v e CC scatters, thereby constraining prediction, 
due to better lateral segmentation and transverse detection. 
Also, improved beam design to further reduce v e from K° . 


Interaction vertex position 


0.00030 


Better lateral segmentation. 


Shower length model 


0.00027 


Better lateral segmentation and transverse detection 
will allow more sophisticated shower identification model. 


Counter efficiency and noise 


0.00023 


Segmented scintillator strips of the type 






developed by MINOS 23 will improve this. 


Energy Measurement 


0.00018 


Better lateral segmentation. 


Charm production, strange sea 


0.00047 


In-situ measurement (TJ [5] . 


Rl 


0.00032 


In-situ measurement [I] [5] . 


rr v Irr" 
G I O 


u.uuuzz 


Likely to be at a similar level. 


Higher Twist 


0.00014 


Recent results reduce this error [24] . 


Radiative Corrections 


0.00011 


New analysis underway, see text below. 


Charm Sea 


0.00010 


Measured in-situ using wrong-sign muon production in DIS. 


Non-isoscalar target 


0.00005 


Glass is isoscalar 



TABLE IV: Source and value of NuTeV errors on sin 2 8w, and reason why the error will be reduced in the PW-style analysis of 
NuSOnG. This paper assumes NuSOnG will reduce the total NuTeV error by a factor of two. This is achieved largerly through 
the improved statistical precision and requires only a 90% reduction in the overal NuTeV systematic error. This table argues 
that a better than 90% reduction is likely, but further study, once the detector design is complete, is required. 



ticlcs which are complementary to the APV and Miller- 
scattering corrections. Thus, these results may not be 
in conflict with NuTeV. The NuSOnG measurement will 
provide valuable additional information on this question. 

Since the NuTeV result was published, more than 
300 papers have been written which cite this result. 
Several "Standard-Model" explanations have been sug- 
gested. While some constraints on these ideas can come 
from outside experiments, it will be necessary for any fu- 
ture neutrino scattering experiment, such as NuSOnG, 
to be able to directly address these proposed solutions. 
Also various Beyond Standard Model explanations have 
been put forward; those which best explain the result 
require a follow-up experiment which probes the neu- 
tral weak couplings specifically with neutrinos, such as 
NuSOnG. Here, we consider the explanations which are 
"within the Standard Model" and address the Beyond 
Standard Model later. 

Several systematic adjustments to the NuTeV result 
have been identified since the result was published but 
have not yet been incorporated into a new NuTeV anal- 
ysis. As discussed here, the corrections due to the two 
new inputs, a new K e3 branching ratio and a new strange 
sea symmetry, are significant in size but are in opposite 
direction - away and toward the Standard Model. So a 
re-analysis can be expected to yield a central value for 
NuTeV which will not change significantly. However, the 
error is expected to become larger. 

In 2003, a new result from BNL865 [29] yielded a K e3 
branching ratio which was 2.3ct larger than past mea- 
surements and a value of \V US \ 2 which brought the CKM 
matrix measurements into agreement with unitarity in 
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FIG. 7: Effect of various "Standard Model" explanations on 
the NuTeV anomaly. The t/-axis is the deviation (S sin 2 6w = 
sin 2 6»^ M - sin 2 0$" TeV ). The solid line is the published 
NuTeV deviation. Thick black lines extending from the 
NuTeV deviation show the range of possible pulls from NLO 
QCD and various isospin violation models. Note that the 
isospin violation models are mutually exclusive and so should 
not be added in quadrature. They are, from left to right, the 
full bag model, the meson cloud model, and the isospin QED 
model. 
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the first row [3U]. The measurement was confirmed by 
CERN NA48/2 [31]. The resulting increased K e3 branch- 
ing ratio [12] increases the absolute prediction of intrin- 
sic v e s in the NuTeV beam. This does not significantly 
change the error because the error on Ke3 was already 
included in the analysis. However, it introduces a cor- 
rection moving the NuTeV result further away from the 
Standard Model, since it implies that in the original anal- 
ysis, NuTeV under-subtracted the v e background in the 
NC sample. The shift in sin 2 8w can be estimated in a 
back of envelope calculation to be about ^0.001 away 
from the Standard Model [32"] . 

The final NuTeV measurement of the difference be- 
tween the strange and anti-strange sea momentum dis- 
tributions, was published in 2007 [33 . This "strange sea 
asymmetry" is defined as 



xs (x) = xs(x) — xs(x), 



(33) 



Because of mass suppression for the production of charm 
in CC scatters from strange quarks, a difference in 
the momentum distributions will result in a difference 
in the CC cross sections for neutrinos and antineutri- 
nos. Thus a correction to the denominator of Eq. (25) 



would be required. The most recent next-to-leading or- 
der analysis finds the asymmetry, integrated over x is 
0.00195 ±0.00055 ±0.00138 [33,. An integrated asymme- 
try of 0.007 is required to explain the published NuTeV 
result |33j . and so one can estimate that this is a shift 
of about 0.0014 in sin 2 9w toward the Standard Model. 
In this case, the errors on the NuTeV result will become 
larger because this effect was not originally considered 
in the analysis. A very naive estimate of the size of the 
increase can be derived by scaling the error on the in- 
tegrated strange sea, quoted above, and is about 0.001 
toward the Standard Model. If this naive estimate of 
the systematic error is borne out, then this could raise 
the NuTeV error on sin 2 9 W from 0.0016 to 0.0018. Nu- 
SOnG will directly address the strange sea asymmetry in 
its QCD measurement program, as described in ref. [2J. 

In ref. [34], additional electromagnetic radiative cor- 
rections have been suggested as a source of the discrep- 
ancy. However, this paper only considered the effect of 
these corrections on R v and not R v and for fixed beam 
energy of E v = 80 GeV. The structure of the code from 
these authors has also made it difficult to modify for use 
in NuTeV. This has prompted a new set of calculations 
by other authors which are now under way |16j . There 
are, as yet, only estimates for the approximate size of 
newly identified effects, which are small. 

The NuTeV analysis was not performed at a full NLO 
level in QCD; any new experiment, such as NuSOnG will 
need to undertake a full NLO analysis. This is possible 
given recently published calculations [35] [35], including 
those on target mass corrections [37]. On Fig. 7] we show 
an early estimate of the expected size and direction of 
the pull [38]. On this plot, the solid horizontal line indi- 
cates the deviation of NuTeV from the Standard Model. 
The thick vertical lines, which emanate from the NuTeV 



deviation, show the range of pulls estimated for various 
explanations. The range of pull for the NLO calculation 
is shown on the left. 

The last possibility is that there is large isospin vi- 
olation (or charge symmetry violation) in the nucleus. 
The NuTeV analysis assumed isospin symmetry, that is, 
u(x) p = d(x) n and d(x) p = u(x) n . Isospin violation can 
come about from a variety of sources and is an interest- 
ing physics question in its own right. NuSOnG 's direct 
constraints on isospin violation are discussed in ref. [2], 
which also considers the constraints from other exper- 
iments. Various models for isospin violation have been 
studied and their pulls range from less than la away from 
the Standard Model to ~ ler toward the Standard Model 
[39"] . We have chosen three examples [39 for illustration 
on Fig. [7j the full bag model, the meson cloud model, 
and the isospin QED model. These are mutually exclu- 
sive models, so only one of these can affect the NuTeV 
anomaly. 



IV. THE TERASCALE PHYSICS REACH OF 
NUSONG 



Even when new states are too heavy to be produced 
at resonance in collisions they can make their presence 
known indirectly, as virtual particles which affect SM 
processes through interference with SM contributions to 
amplitudes. The new heavy states induce small shifts 
in observables from SM predictions, and conversely pre- 
cise measurements of these observables can constrain or 
detect new physics at mass scales well above the ener- 
gies of the colliding particles. In this way the preci- 
sion neutrino scattering measurements at NuSOnG will 
place TeV-scale indirect constraints on many classes of 
new physics, or perhaps detect new physics by measur- 
ing deviations from SM predictions. The effects of new 
high-scale physics may be reduced to a small number of 
effective operators along with corresponding parameters 
which may be fit to data. Although the particular set 
of operators used depends on broad assumptions about 
the new physics, the approach gives a parameterization 
of new physics which is largely model-independent. 

For concrctencss we will assume that NuSOnG will be 
able to measure the neutrino ES/IMD ratio to a preci- 
sion of 0.7%, a(D^e) (normalized as per Sec. Ill B I to 
1.3%, and that NuSOnG will be able to halve the errors 
on NuTeV's measurement of DIS effective couplings, to 
Agl = 0.0007 and Ag 2 R = 0.0006 (where g L and g R were 



defined in Eqs. (29 1 and (30)) 



We first parameterize new physics using the oblique 
parameters ST, which is appropriate when the impor- 
tant effects of the new physics appear in vacuum polar- 
izations of gauge bosons. We next assume new physics 
effects manifest as higher-dimensional operators made of 
SM fermion fields. We separately consider the possibil- 
ity that the gauge couplings to neutrinos are modified. 
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Topic 


Contribution of NuSOnG Measurement 


Oblique Corrections 


Four distinct and complementary probes of S and T. 

In the case of agreement with LEP/SLD: ~25% improvement in electroweak precision. 


Neutrino-lepton NSIs 


Order of magnitude improvement in neutrino-electron effective couplings measurements. 
Energy scale sensitivity up to ~ 5 TeV at 95% CL. 


Neutrino-quark NSIs 


Factor of two improvement in neutrino-quark effective coupling measurements. 
Energy scale sensitivity up to ~ 7 TeV at 95% CL. 


Mixing with Neutrissimos 


30% improvement on the e-family coupling in a global fit. 
75% improvement on the ^-family coupling in a global fit. 


Right-handed Couplings 


Complementary sensitivity to gn/gL compared to LEP. 

Order of magnitude improvement compared to past experiments. 



TABLE V: Summary of NuSOnG's contribution to general Terascale physics studies. 



Realistic models usually introduce several new operators 
with relations among the coefficients; we consider several 
examples. A summary of the contributions of NuSOnG 
to the study of Terascale Physics is provided in Table IV] 



A. Oblique corrections 

For models of new physics in which the dominant loop 
corrections are vacuum polarization corrections to the 
SU{2) l x U(1)y gauge boson propagators ("oblique" cor- 
rections), the STU [301 SI] parameterization provides a 
convenient framework in which to describe the effects of 
new physics on precision electroweak data. Differences 
between the predictions of a new physics model and those 
of a reference Standard Model (with a specified Higgs bo- 
son and top quark mass) can be expressed as nonzero val- 
ues of the oblique correction parameters S, T and U. T 
and U are sensitive to new physics that violates isospin, 
while S is sensitive to isospin-conserving physics. Pre- 
dictions of a Standard Model with Higgs or top masses 
different from the reference Standard Model may also be 
subsumed into shifts in S and T (in many models U is 
much smaller than S and T and is largely unaffected by 
the Higgs mass, so it is often omitted in fits). Within a 
specific model of new physics the shift on the ST plot 
away from the SM will be calculable [42,. For example, 

• A heavy Standard Model Higgs boson will make 
a positive contribution to S and a larger negative 
contribution to T. 

• Within the space of Z' models, a shift in almost 
any direction in ST space is possible, with larger 
shifts for smaller Z' masses. 

• Models with a fourth-generation of fermions will 
shift S positive, and will shift T positive if there 
are violations of isospin. 

In constructing models incorporating several types of new 
physics the corresponding shifts to S and T combine; if 
contributions from different sectors are large, then they 
must conspire to cancel. 



EH 0.0 




FIG. 8: The impact of NuSOnG on the limits of S and T. 
The reference SM is m t = 170.9 GeV, and ran = 115 GeV. 
la bands due to NuSOnG observables are shown against the 
90% contour from LEP/SLD. The central ellipses are the 68% 
and 90% confidence limit contours with NuSOnG included. 



See Eqs. (29 1 and (301 for the definitions of gh and gu. 



The constraints on S and T from the full set of preci- 
sion electroweak data strongly restrict the models of new 
physics which are viable. The strongest constraints are 
from LEP/SLD, which give a current bound of 



S = -0.02 ±0.11 
T = +0.06 ±0.13 
Corr(S, T) = 0.91. 



(34) 



The ES and DIS measurements from NuSOnG provide 
four distinct and complementary probes of S and T, as 
shown in Fig. [8] If the target precision is achieved, and 
assuming the NuSOnG agree with SM predictions, Nu- 
SOnG will further reduce the errors on S and T from the 
LEP/SLD values to 



S 
T 



-0.05 ±0.09 
±0.02 ±0.10 
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Corr(S,T) = 0.87 



(35) 



where 



The ~ 25% reduction in the errors is primarily due to the 
improved measurement of g\. We note that the error g\ 
is likely to be further reduced (see Sec. HIE), and so the 
this is conservative estimate of NuSOnG's contribution 
to the physics. 



9v 


= 2gl9 f v = 


9a 


= 2gl9 f A = 


9l 


= WlqI = 


9r 


= 2gl9 f R = 



(iS) . 



p(-Q* sin 2 6 



(38) 



B. Non-standard interactions 

NuSOnG will probe new physics that modifies 
neutrino-quark and neutrino-electron scattering. If the 
masses associated to the new degrees of freedom are much 
larger than the center of mass energy (s = 2m e i?bcam 
0.5 GeV 2 ) then modifications to these processes are well- 
described by higher-dimensional effective operators. In 
the context of neutrino reactions, these operators are also 
referred to as non-standard interactions (NSI's). In a 
model-independent effective Lagrangian approach these 
effective operators are added to the SM effective La- 
grangian with arbitrary coefficients. Expressions for ex- 
perimental observables can be computed using the new 
effective Lagrangian, and the arbitrary coefficients can 
then be constrained by fitting to data. Typically, bounds 
on the magnitude of the coefficients are obtained using 
only one or a few of the available effective operators. This 
approach simplifies the analysis and gives an indication of 
the scale of constraints, although we must be mindful of 
relationships among different operators that will be im- 
posed by specific assumptions regarding the underlying 
physics. 

To assess the sensitivity of NuSOnG to "heavy" new 
physics in neutral current processes, we introduce the 
following effective Lagrangian for neutrino-fermion inter- 
actions m Ha m : 



c 



NSI — 



-V2G F 
-2y/2G F 



fA 



(36) 



where a, (3 — e,[i,T and L, R represent left-chiral and 
right-chiral fcrmion fields. If a ^ /?, then the a /3 
terms must be Hermitian conjugates of each other, i.e. 



£(3o 



'a/3 



NuSOnG is sensitive to the /? = /i couplings. 



This effective Lagrangian is appropriate for parameteriz- 
ing corrections to neutral current processes; an analysis 
of corrections to charged-current processes requires a dif- 
ferent set of four-fermion operators. 

Assuming e a p = for a ^ /3 we need consider only the 
terms e£* (* = V, A, L, R). If we rewrite Eq. ^ as 



C 



-V2G F 
-2a/2G f 



v^^Plv 

Vlp,PhV 



stf f-ff-g"! /V75/ 

gl f frp L f 

+ g R f fi fl P R f 



(37) 



then we see that adding Eq. ([36j to the SM Lagrangian 
will simply shift the effective couplings: 



9v 


— ► g v 


= g v - 


he 


g/ 


~vf 

— ► 9 A 


= g/ - 


he 


g L 


~vf 

— ► 9l 


= g L - 


he 


g R 


~vf 

— ► 9r 


= g R - 


he 



fv 

fA 
MM 
fL 
Hfi ■ 
fR 
mm ' 



(39) 



Consequently, errors on the g P s translate directly into 
errors on the e^'s, P = V, A or P = L, R. 



1. Neutrmo-lepton NSI 

A useful review of present constraints on non-standard 
neutrino-electron interactions can be found in ref. [45] . 
As this paper states, and as we show below, an improved 
measurement of neutrino-elecron scattering is needed. 

The world average value for neutrino-electron effective 
couplings, dominated by CHARM II, is 



9v 

g7 



Corr( 5 -, 5 ^) = 



-0.040 ±0.015 
-0.507 ±0.014 
-0.05 . 



(40) 



The current la bounds from CHARM II, Eq. (40 1 trans- 
lates to \e e J^\ < 0.01, (P — L.R) with a correlation 
of 0.07 jH]. At the current precision goals, NuSOnG's 
v^e and V^e will significantly reduce the uncertainties on 
these NSI's, to 



'MM ' MM ' 



-eV\ 
'MM I 
-eA\ 
'MM I 
eA\ 



< 0.0036 , 

< 0.0019 , 
= -0.57, 



Corr(e 

or in terms of the chiral couplings, 

|e^| < 0.0015 
|e^| < 0.0025 
Corr(e^,e^) = 0.64.. 



(41) 



(42) 



Even in the absence of a a{y^e) measurement e^h and 
ef$ can be constrained from the v^e scattering data alone 
through a fit to the recoil electron energy spectrum (see 
Eq. (§}). 

We first consider the constraint on ef,^ and ef/? from 
the total cross section ayv^e). It is convenient to recast 
the effective interaction slightly, as 
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-NSI 



-2V2G F 
V2 



al.PL^ e±&fP L e + e^erfP R e 



A 2 



cos 9 ej° Pl & + sin 9 e-f" P^e 
I 



(43) 



The new physics is parameterized by two coefficients A 
and 9. A represents the broadly-defined new physics scale 
while 9 E [0, 2ir] defines the relative coupling of left-chiral 
and right-chiral electrons to the new physics. As an ex- 
ample, a scenario with a purely "left-handed" Z' that 
couples to leptons with coupling g' would be described 
by A oc Mz'/g' and 9 = or 9 = n, depending on the 
relative sign between g' and the clcctrowcak couplings. 
A and 9 are related to to the NSI parameters in Eq. ( 36 ) 

by 



cos f 



2G F A 2 ' 



eR _ _ 



sin ( 



2G F A 2 ' 



(44) 



Note that we have generalized from our assumption of 
the previous section and not taken a — fi necessarily. At 
NuSOnG, new physics modifies (pseudo)elastic neutrino- 
electron scattering. Here we use the word "pseudo" to 
refer to the fact that we cannot identify the flavor of the 
final-state neutrino, which could be different from the 
incoming neutrino flavor in the case of flavor changing 
neutral currents. 

The shift in the total cross section is 

M"/*<0 { 2 9l eft + «0 2 } + H 2 9 V R e tffr + (^') 2 } 



aiv^e) 



o?r) 2 + §(5«) 5 



516 GeV 
A 



+0.096 ( 516 ^ V ) (l + 2cos 2 , 



(45) 



where 



tan ( 



9j[ 



-0.28 



(46) 



When 0(e 2 ) terms are negligible, a 0.7% measurement 
of a(iy^e) translates into a 95% confidence level bound of 



A > 



(4.4 TeV) x a/|cos(0-0)| 



(47) 



from elastic scattering. 

The measurement of the electron recoil energy will al- 
low us to do better. Fig. [9jdark line) depicts the 95% 
confidence level sensitivity of NuSOnG to the physics 
described by Eq. (43 1 when v a = obtained af- 
ter fitting the recoil electron kinetic energy distribu- 
tion. Fig. [^closed contour) represents how well NuSOnG 
should be able to measure A and 9, at the 95% level. 
Weaker bounds from pseudoelastic scattering are also 




100 150 200 

6 (degrees) 



FIG. 9: (DARK LINES) 95% confidence level sensitivity of 
NuSOnG to new heavy physics described by Eq. ( 43 1 when 



fa = (higher curve) and u a / z/ M (lower curve). (CLOSED 
CONTOURS) NuSOnG measurement of A and 0, at the 95% 
level, assuming v a — i/ M , A = 3.5 TeV and 9 = 27r/3 (higher, 
solid contour) and v a 7^ z/ M , A = 1 TeV and 8 = 4n/3 (lower, 
dashed contour). Note that in the pseudoelastic scattering 
case [y a 7^ v^) 8 and n + 8 are physically indistinguishable. 



scattering. While there will be fewer of these events, they 
should qualitatively improve our ability to pin down the 
new physics parameters given the distinct dependency on 
g^ e and g" A e (see Sec. [Big) . 

Eq. ( 43 1 does not include all effective dimension-six 



operators that contribute to neutrino-electron (pseudo) 
elastic scattering. All neglected terms will either not con- 
tribute at NuSOnG, or were assumed to be suppressed 
with respect to Eq. (43). In turn, terms proportional 



to a right-handed neutrino current PrJu^r lead to neg- 
ligibly small effects since neutrino masses are negligibly 
small and we are dealing with neutrino beams produced 
by pion and muon decay (i.e., for all practical purposes, 
we have a purely left-handed muon neutrino beam and a 
purely right-handed muon antineutrino beam) . Chirality 
violating effective operators (e.g. {vRVL,){eLe-R)) , on the 
other hand, are expected to be suppressed with respect 
to Eq. ( 43 ) by terms proportional to neutrino masses and 



the electron mass (measured in units of A). The reason 
is that, in the limit of masslcss neutrinos or a massless 
electron, chiral symmetry is restored while such operators 
explicitly violate it. For the same reason, dimension- five 
magnetic moment-type operators {yijpavF 1 " 7 ) have also 
been neglected. 



shown. We have not included 



'data" from ^-electron 



We note also that Eq. (43 1 violates SU(2)l unless 
one also includes similar terms where «-> £^ (I — 
e, jJL, t). In this case, certain flavor combinations would be 
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severely constrained by electron-electron scattering and 
rare muon and tau decays. One way around such con- 



straints is to postulate that the operators in Eq. ( 43 ) are 
dimension-eight operators proportional to LH*"f a LH , 
where L is the left-chiral lepton doublet and H is the 
Higgs scalar doublet. In this case, 1/A 2 should be re- 
placed by v 2 /A 4 , where v = 246 GeV is the scale of 
electroweak symmetry breaking. 

Finally, another concern is whether modifications to 
the charged current neutrino-electron (pseudo)quasi- 
elastic scattering ((pseudo)IMD, v^e — * v a \\i) can render 
the translation of NuSOnG data into constraints or mea- 
surements of 6 and A less straightforward. This turns 
out not to be the case, since new physics contributions 
to v^e — ► v a n are already very well constrained by preci- 
sion studies of muon decay. Hence, given the provisos of 
the two previous paragraph, Eq. (43 1 is expected to cap- 



ture all "heavy" new physics effects in (pseudo)elastic 
neutrino electron scattering. 



2. Neutrino- quark NSI 

We next consider the f = u, d case. The change in the 
parameters g\ and g\ (see Eqs. ( 29|30 |) due to the NSI's 



important, as is the case in some of the scenarios we treat 
later, a new analysis is necessary and the bounds above 
cannot be used. This is to be contrasted to the neutrino- 
lepton case, discussed in the previous subsection. 



Neutrissimos, Neutrino Mixing and Gauge 
Couplings 



0.02 




*tfk 



n vu uL i n vd dL 

z 9l e / i A t + z 9l £ fitx 
+0.69 e%- 0.85 e% 



-0.31 e% + 0.15 e a ^ 



-dR 



(48) 



so only these linear combinations are constrained. The 
bounds from NuTeV (rescaled to la bounds from 
ref. jH]) are: 



e«u = -0.0053 ± 0.0020 , 

= +0.0043 ± 0.0016 , 

\eH\ < 0.0035 , 

\e%\ < 0.0073 . 



(49) 



These bounds are obtained by setting only one of the 
parameters be non-zero at a time. If NuSOnG reduces 
the errors on the NuTeV measurement of g\ and g\ by 
a factor of 2, the la bounds on the NSI parameters are 
similarly reduced: 



-fj.fi, I 

-llfll 
-flfl I 



0.001 , 
0.0008 , 
0.002 , 
0.004 . 



(50) 



In terms of a new physics scale defined as A = 1/^/2 Gp£, 
these constraints range from A > 3TeVtoA > 7TeV. 

We note that neutrino-quark scattering will also be 
sensitive to NSIs which correct CC interactions. These 
interactions are not included in Eq. ( 36 1 . If they are 



FIG. 10: Potential constraint on e e and e M from NuSOnG 
(see Eq. (55 l). This is a two-dimensional projection of a 4 



90% CL contour of a fit to all the charge current particle decay 
data + LEP/SLD. 

In those classes of models which include moderately 
heavy electroweak gauge singlet ( "neutrissimo" ) states, 
with masses above 45 GeV, the mixing of the SU (2) in- 
active neutrinos and the sterile states may lead to a sup- 
pression of the neutrino-gauge couplings. The resulting 
pattern of modified interactions is distinct from those 
of the previous section since they will also induce corre- 
lated shifts to the charged-current coupling. For exam- 
ple, Ref. [IB] presents models with one sterile state per ac- 
tive neutrino flavor and intergenerational mixing among 
neutrinos. In these models the flavor eigenstates are lin- 
ear combinations of mass eigenstates, and those mass 
eigenstates too heavy to be produced in final states re- 
sult in an effective suppression of the neutrino-gauge bo- 
son coupling. This suppression may be flavor-dependent 
depending on the structure of the neutrino mixing ma- 
trix. If the mass matrix contains Majorana terms, such 
models permit both lepton flavor violation and lepton 
universality violation. 

Neutrinos couple to the W and the Z through interac- 
tions described by: 



-°-w- 

V2 " 
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zsc 



(51) 



where t — e, /i,r. If the neutrinos mix with gauge sin- 
glet states so that the SU(2)l interaction eigenstate is a 
superposition of mass eigenstates fought and i^heavy 

vil = vi, light cos 9 e + ^,hcavy sin 9 e , (52) 

then the interaction of the light states is given by 

C = 

( -?-W- trfviMtt + vtp&rflt J cos 5/ 



+ (^J^ ^,Ught7 M ^,lightJ COS 2 0i . 
Defining 

Q = 1 — cos 2 Of ■ 
the shift in the Lagrangian due to this mixing is 



(53) 



(54) 



Q 
2 



e 
2^ 



(55) 



where we have dropped the subscript "light" from the 
neutrino fields. 

Lepton universality data from W decays and from 
charged current tt,t and K decays [37] constraint differ- 
ences eg i — q. . LEP/SLD and other precision electroweak 
data will imposed additional constraints on q in combi- 
nation with the oblique parameters, as will NuSOnG. A 
fit to all the charge current decay data and LEP/SLD 
with S, T, e e and e„ yields 



S = -0.05 ±0.11 , 

T = -0.44 ±0.28 , 

e e = 0.0049 ± 0.0022 , 

e u = 0.0023 ± 0.0021 . 



(56) 



If we now included hypothetical data from NuSOnG, 
assuming NuSOnG achieves its precision goals and mea- 
sures central values consistent with the Standard Model, 
we see the constraints on e M and e e are substantially im- 
proved. In this case, the fit yields 



S = 0.00 ±0.10 , 

T = -0.11 ±0.12 , 

e e = 0.0030 ±0.0017 

e„ = 0.0001 ±0.0012. 



fx — > g^/ implies e e e M ss 0. This type of model has been 
proposed as a solution to the NuTeV anomaly. If we 
take take only one of e e or to be nonzero (to respect 
the constraint from \i — * cy), the NuTeV value of g\ is 
accommodated in the fit by best-fit values of e that are 
large and positive and best-fit values of T are large and 
negative (consistent with a heavy Higgs) . 



D. Right-handed coupling of the neutrino to the Z 

In the Standard Model, neutrino couplings to the W- 
and Z-bosons are purely left-handed. The fact that the 
neutrino coupling to the Ty-boson and an electron is 
purely left-handed is, experimentally, a well-established 
fact (evidence includes precision measurements of pion 
and muon decay, nuclear processes, etc.). By contrast, 
the nature of the neutrino coupling to the Z boson is, 
experimentally, far from being precisely established [SU] . 
The possibility of a right-handed neutrino- Z-boson cou- 
pling is not included in the previous discussions, and is 
pursued separately in this subsection. 

The best measurement of the neutrino coupling to the 
Z-boson is provided by indirect measurements of the in- 
visible Z-boson width at LEP. In units where the Stan- 
dard Model neutrino- Z-boson couplings are g v L = 0.5, 
cf R = 0, the LEP measurement [FT] translates into 
(9l) 2 + (9r) 2 = °- 2487 ± 0.0010. Note that this result 
places no meaningful bound on g v R . 

Precise, model-independent information on g u L can 
be obtained by combining + e scattering data from 
CHARM II and LEP and SLD data. Assuming model- 
independent couplings of the fermions to the Z-boson, 

+ e scattering measures g v L = yfp/2, while LEP and 
SLD measure the left and right-handed couplings of the 
electron to the Z. The CHARM II result translates into 
\g v L \ = 0.502 ± 0.017 [5U], assuming that the charged- 
current weak interactions produce only left-handed neu- 
trinos. In spite of the good precision of the CHARM II 
result (around 3.5%), a combination of all available data 
allows \g V ]il g v j\ ~ 0.4 at the two a confidence level [5U] , 

Significant improvement in our understanding of g v R 
can only be obtained with more precise measurements of 
v+e scattering, or with the advent of a new high intensity 
e + e~ collider, such as the ILC. By combining ILC run- 
ning at the Z-boson pole mass and at y/s = 170 GeV, 
\9r/9l\ ~ 0.3 could be constrained at the two a level 



(57) 



after analyzing e + e~ — > 7+missing energy events |50j . 
Assuming that g v L can be measured with 0.7% uncer- 
Fig. |10| shows the two dimensional cross section in the tainty, Fig. 11 depicts an estimate of how precisely g v R 

could be constrained once NuSOnG "data" is combined 
Fig. [TT|left) considers the hypothesis 
In 



e e -e M plane of the four dimensional fit. The likelihood 



contours are 2D projections. Though not obvious from 
the figure, it is NuSOnG 's improved measurement of g\ 
which contributes the most to strengthening the bounds 
on the €£. 

In models of this class lepton flavor violating decays 
such as — ► cy impose additional constraints on prod- 
ucts e^ef. . For example, the strong constraint from 



with LEP data. 

that the Standard Model expectations are correct, 
this case, NuSOnG data would reveal that gn/gh is less 
than 0.2 at the two sigma level. On the other hand, 
if gnlgh = 0.25 - in good agreement with the current 
CHARM II and LEP data - NuSOnG data should reveal 
that gji 7^ at more than the two sigma level, as depicted 
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FIG. 11: Precision with which the right-handed neutrino-Z-boson coupling can be determined by combining NuSOnG mea- 
surements of g v L with the indirect determination of the invisible Z-boson width at LEP if (left) the v + e scattering measurement 
is consistent with the Standard Model prediction g v h = 0.5 and (right) the v + e scattering measurement is significantly lower, 
g v L = 0.485, but still in agreement with the CHARM II measurement (at the one sigma level). Contours (black, red) are one 
and two sigma, respectively. The star indicates the Standard Model expectation. 



in Fig. 11 Vight) 



The capability of performing this measurement in 
other experiments has been examined. The NuSOnG 
measurement compares favorably, and complements, the 
ILC capabilities estimated in [50] . Ref [52] studied mea- 
surements using other neutrino beams, including reactor 
fluxes and beta beams. NuSOnG's reach is equivalent to 
or exceeds the most optimistic estimates for these various 
neutrino sources. 



V. SPECIFIC THEORETICAL MODELS AND 
EXPERIMENTAL SCENARIOS 

If NuSOnG's measurements agree with the SM within 
errors, we will place stringent constraints on new physics 
models; if they disagree, it will be a signal for new 
physics. In the latter case the availability of both DIS 
and ES channels will improve our ability to discriminate 
among new physics candidates. NuSOnG will also pro- 
vide an important complement to the LHC. The LHC 
will provide detailed information about the spectrum of 
new states directly produced. However, measurements of 
the widths of these new states will provide only limited 
information about their couplings. NuSOnG will probe 
in multiple ways the couplings of these new states to neu- 
trinos and to other SM particles. 

In this section we provide several case studies of 
NuSOnG sensitivity to specific models of new physics. 
These include several typical Z' models, leptoquark mod- 
els, models of R-parity violating supersymmetry, and 
models with extended Higgs sectors. We examine how 



these will affect v^e ES and u^N DIS at tree-level. Our 
list is far from exhaustive but serves to illustrate the pos- 
sibilities. We summarize our contributions in Table IVl 



The opposite way to approach this problem is to ask: 
in the face of evidence for new Terascale Physics, how 
can we differentiate between specific models? NuSOnG 
has the potential to discover new physics through indi- 
rect probes, in the event that one or more of its mea- 
surements definitively contradicts SM predictions. We 
discuss several possible patterns of deviation of model- 
independent parameters from SM predictions and some 
interpretations in terms of particular models. This is pre- 
sented in the context of various expectations for LHC to 
illustrate how NuSOnG enhances the overall physics pro- 
gram. Since the NuTeV reanalysis is ongoing, and since 
the ES constraints from CHARM-II are weak, it is pru- 
dent that we commit to no strong assumptions about the 
central value of the NuSOnG measurements but instead 
consider all reasonable outcomes. 



A. Sensitivity in the Case of Specific Theoretical 
Models 

We next consider the constraints imposed by the pro- 
posed NuSOnG measurements on explicit models of BSM 
physics. An explicit model provides relations among 
effective operators which give stronger and sometimes 
better-motivated constraints on new physics than is ob- 
tained from bounds obtained by considering effective op- 
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Model 


Contribution of NuSOnG Measurement 


Typical Z' Choices: (B — xL),(q — xu),(d + xu) 


At the level of, and complementary to, LEP II bounds. 


Extended Higgs Sector 


At the level of, and complementary to r decay bounds. 


R-parity Violating SUSY 


Sensitivity to masses ~ 2 TeV at 95% CL. 
Improves bounds on slepton couplings by ~ 30% and 
on some squark couplings by factors of 3-5. 


Intergenerational Leptoquarks with non-degenerate masses 


Accesses unique combinations of couplings. 

Also accesses coupling combinations explored by n decay bounds, 
at a similar level. 



TABLE VI: Summary of NuSOnG's contribution in the case of specific models 



„ 6 
> 

w 5 



Sensivity Range for NuSOnG at >95% CL 





U(1) B - x l 






u(i)d~xu 


en 


—x 
—x 


-1 

-(2 + sc)/3 


x/3 
-1/3 


(-l + x)/3 
x/3 



TABLE VII: Charges of u^l ,eL,en under 4 phenomenolog- 
ically viable classes of U(l)' symmetries. Each value of x 
corresponds to a different U(l)' symmetry that is considered. 



FIG. 12: Some examples of NuSOnG's 2<r sensitivity to new 
high-mass particles commonly considered in the literature. 
For explanation of these ranges, and further examples, see 
text. 



erators one by one, but at the expense of the generality 
of the conclusions. Many models can be analyzed using 
the effective Lagrangian of Eq. ( 36 1 , but others introduce 



new operators and must be treated individually. The list 
of models considered is not exhaustive, but rather illus- 
trates the new physics reach of NuSOnG. 



1. Z' models 

Massive Z' fields are one of the simplest signatures of 
physics beyond the Standard Model. (For a recent re- 
view, see [53 .) Z' vector bosons are generic in grand 
unified theories and prevalent in theories that address 
the electroweak gauge hierarchy. They may stabilize 
the weak scale directly by canceling off quadratic diver- 
gences of Standard Model fields, as in theories of extra- 
dimensions or Little Higgs theories. In supersymmetric 
models, Z' fields are not needed to cancel quadratic di- 
vergences, but are still often tied to the scale of soft- 
breaking (and hence the electroweak scale). In these last 
two cases, the Z' typically has a TeV-scale mass, and is 
an attractive target for NuSOnG. 



If the Z' mass is sufficiently large, its exchange is well- 
described at NuSOnG energies by the effective operator 
of Eq. (43 1. In this case, the new physics scale is related 
to the Z^model by A ~ M Z '/gz', the ratio of the Z' 
mass to its gauge-coupling. Further model-dependence 
shows up in the ratio of fermion charges under the U(l)' 
symmetry associated with the Z' , and the presence of 
any Z — Z' mixing. With reasonable theoretical assump- 
tions, the absence of new sources of large flavor-changing 
neutral currents, the consistency of Yukawa interactions, 
and anomaly cancellation with a minimal number of ex- 
otic fermions, the number of interesting models can be 
reduced substantially, to four discrete families of generic 
(7(1)' models each containing one free parameter, x [54] . 
In Table V A 1 we indicate the charges of v^hi e L,^R un- 
der these families of U(l)' symmetries. 

Using the sensitivity of NuSOnG to the scale A in 
scattering shown in Figure [9] we can bound the combi- 
nation Mz 1 / gz' for the four families of Z' models as a 
function of a;. It is important to note that these bounds 
are competitive with the LEP-II bounds found in |54j . 
which are based on Z' decays to all fermions, not just 
electrons and neutrinos. 

There are Z' models which distinguish among genera- 
tions can affect neutrino scattering. These will be probed 
by NuSOnG at the TeV scale [S3 S EH EH EH] ■ Among 
these, B — 3L M was suggested as a possible explanation 
for the NuTeV anomaly [SHIED], however, we show here 
that this is not the case. Nevertheless, it remains an 
interesting example to consider. 

In the gauged B - 3L M the Z' modifies v^N DIS. The 
exchange of the Z' between the and the quarks induces 
operators with coefficients 



,uR 



-dR 
'MM 
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Eq. ( 36 ) which with coefficient 



FIG. 13: 95% confidence level sensitivity of NuSOnG to the 
indicated Z' models. The charges of the electrons and neutri- 
nos under the underlying (7(1)' gauge symmetry are described 
in Table |V A 1| The bounds are plotted as functions of the 
parameter x, which scans over allowed fermion charges for 
each family of (7(1)' symmetries, versus the ratio M z //g Z '- 



e eL = - 

MM 



1 M 

V2G F Ml 



e eR = 
mm u 



From Eq. Q, the 95% bound is: 



I A. 



> 5.2 TeV, 



(62) 



(63) 



Cfl 



competitive with current bound from r-decay of 5.4 TeV. 



3. R-parity violating SUSY 

Assuming the particle content of the Minimal Super- 
symmetric Standard Model (MSSM) , the most general Im- 
parity violating superpotential (involving only tri-linear 
couplings) has the form 64J 

where .£,-, Qi, Di, and Ui are the left-handed MSSM 
superfields defined in the usual fashion, and the sub- 
scripts i,j,k= 1,2,3 are the generation indices. SU (2) j, 
gauge invariance requires the couplings Xijk to be anti- 
symmetric in the first two indices: 



2y/2G F Mf 



9% F 

= £B-3L U 



which shift g\ and g\ by 



2s 2 



£B-3L U 



(58) 



(59) 



It should be noted that since Eb-zl^ is negative, this 
shows that both g 2 L and g\ will be shifted positive. This, 
in fact, excludes gauged B — 3L M as an explanation of 
the NuTeV anomaly. With this said, a NuSOnG mea- 
surement of g\ and g 2 R that improves on NuTeV errors 
by a factor of 2 yields a 2er bound 



A 



ijk 



= -A 



(65) 



The purely baryonic operator UiDjDk is irrelevant to 
neutrino scattering, so only the 9 Xijk and 27 A'^ fe cou- 
plings are of interest. 



From the LLE part of the Eq. ( 64 1 slepton exchange 



will contribute to v^e ES at NuSOnG. These induce four- 



fcrmion operators appearing in Eq. ( 36 1 with correspond- 
ing coefficients 



.eh 



'MM 



1 



E 



|A. 



21fc 



1 



M? 

k = l e kR 



|A 2j i| 2 
M? 

3 e it 



(66) 



9z- 



> 2.2 TeV 



(60) 



which is comparable and complementary to the existing 
bound from DO, and thus interesting to consider. 



2. Models with extended Higgs sectors 



If we place bounds on the sleptons one at a time, then 



Eq. (42 1 translates to the 2cr bounds shown in Table VIII 



presented for masses of 100 GeV. To rescale to different 
masses, use ( ioogcV ) ' This can be compared to current 
bounds Ref. [55]. NuSOnG improves all of these bounds. 
From the LQD part of Eq. ( 64 ) , squark exchange will 



contribute to contribute to NC v^N DIS and CC v^N 
DIS. The resulting shifts in g\ and g 2 R are 



In the Zee and Babu-Zee [53] models, an isosinglct 
scalar h + with hypcrcharge Y = +1 is introduced, which 
couples to left-handed lepton doublets I as 



Sg 



Sg 2 R = 2 



ud dR _ 2 

9r e nfi 9r £ c 



(67) 



C-h = X ab (e c aL ia 2 £ bL ) h + + h.c. , 



(61) 



where (ab) are flavor indices: a, b = e, fx, r. The exchange 
of a charged Higgs induces the effective operator from 



where 



3 



^2ifel 2 



dL 

r fe=l d kR 
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Coupling 


yo/o iNUoVjnljr Douna 


current 95% bound 


|Al21 




0.03 


0.05 (Vud) 






0.04 


0.05 (V u d) 


|Al23 




0.04 


0.05 (Vud) 


| A231 




0.05 


0.07 (r decay) 


1^211 




0.05 


0.06 (tv decay) 


1^212 




0.06 


0.06 (n decay) 


1^213 




0.06 


0.06 (-7T decay) 


1^221 




0.07 


0.21 (D meson decay) 


^231 




0.07 


0.45 (Z-> n + n~) 



TABLE VIII: Potential bounds on the R-parity violating LLE 
(top) and LQD (bottom) couplings from NuSOnG, assuming 
that only one coupling is non-zero at a time for each set. All 
squark and slepton masses are set to 100 GeV. To obtain lim- 
its for different masses, rescale by ( kxhScV ) • Current bounds 
are from Ref. 1651. 
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(68) 



e„„ and are associated with terms of Eq. ( 36 1 , while 
e r is associated with a four-fcrmion interaction that cor- 



rects charged currents, 



- 2V2G F e 
The shifts in g\ and g\ are: 

Sgl = 2(gl d 



h.c. 



Sgl 



2 \ 

2 „dL 1 c^vd^dR 



(69) 



(70) 



Assuming the projected precision goals for NuSOnG on 
g\ and g R , and allowing only one of the couplings to be 
nonozero at a time, the 2c bounds are given in Tablc|VTTT] 
mass of 100 GeV, in all cases. To obtain limits for differ- 
ent masses, one simply rescales by ( 100 A g cV ). NuSOnG's 
measurements are competitive with ir decay bounds, and 
improves the current bounds on the 221 and 231 cou- 
plings by factors of 3 and 5, respectively. 



4- Intergenerational leptoquark models 

Measurements of g\ and g R are sensitive to lepto- 
quarks. Because the exchange of a leptoquark can in- 
terfere with both W and Z exchange processes, we can- 
not use the limits on the NSI's of Eq. (36), since we 



must also include the effects of the four-fermion opera- 
tors associated with charged-current processes. Instead, 
the interactions of leptoquarks with ordinary matter can 
be described in a model-independent fashion by an effec- 
tive low-energy Lagrangian as discussed in Refs. [66] HE] 
for generation-universal leptoquark couplings. For lepto- 
quarks to contribute to v^N DIS, they must couple sec- 
ond generation leptons to first generation quarks, so we 



use the more general Lagrangian of jHZl [55] , which allows 
the coupling constants to depend on the generations of 
the quarks and leptons that couple to each leptoquark. 
We summarize the quantum numbers and couplings of 



the various leptoquarks fields in Table IX our notation 
conventions are those of Ref. [69] . 



The four-fermion operators induced by leptoquark ex- 
change will affect NC and/or CC processes, and at Nu- 
SOnG the effect manifests itself in shifts g\ and g R . As- 
suming degenerate masses within each iso-multiplet, the 
shifts in g\ and g R can be written gcncrically as 



Sgl 



Sgl 



= c. 



= c. 



|AFg| 2 /M| Q 

g 2 /M^ 

lA^P/M 2 



\12 12 
A LQ\ 



LQ 



C L 
C R 



LQ I 



3 2 /M 



w 



AV2G F M 2 LQ 



(71) 



where \ l 2 q denotes the (ij) — (12) coupling of the lepto- 
quark and Mlq is its mass. In table [X we list what 
they are, and in figure [14] we plot the ^dependence of 
Sgl and Sg R on the ratio |Alq| 2 /M|q. Table |x| also 
lists the projected NuSOnG bounds on the coupling con- 
stants [70] . Existing bounds on Si, S3, V\, and V3 cou- 
plings from R n — Br(it — > ev)/Br(7r — > [iv) are already 
much stronger, but could be circumvented for S3 and 
V3 if the masses within the multiplct are allowed to be 
non-degenerate. 



B. Interplay with LHC to Isolate the Source of 
New Physics 

By the time NuSOnG runs, the LHC will have accu- 
mulated a wealth of data and will have begun to change 
the particle physics landscape. The message from LHC 
data may be difficult to decipher, however. As discussed 
below, NuSOnG will be able to help elucidate the new 
physics revealed at the LHC. The discovery of a Higgs 
along with the anticipated measurement of the top mass 
to 1 GeV precision would effectively fix the center of the 
ST plot and will enhance the power of the precision elec- 
troweak data as a tool for discovering new physics. If ad- 
ditional resonances are discovered at the LHC, it is still 
likely that little will be learned about their couplings. 

The NuSOnG experiment provides complementary in- 
formation to LHC. Rather than generalize, to illustrate 
the power of NuSOnG, two specific examples are given 
here. We emphasize that these are just two of a wide 
range of examples, but they serve well to demonstrate 
the point. Here we have chosen examples from typical 
new physics models other than Z' models which were dis- 
cussed above, in order to demonstrate the physics range 
which can be probed by NuSOnG. 
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TABLE IX: Quantum numbers of scalar and vector leptoquarks with SU(3)c x SU(2)l X U(1)y invariant couplings to quark- 
lepton pairs (Q cm = I 3 + Y) [T^] ■ 
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TABLE X: Potential and existing 95% bounds on the leptoquark couplings squared when the leptoquark masses are set to 
100 GeV. To obtain the limits for different leptoquark masses, multiply by (Mlq /100 GeV) 2 . Existing bounds on the Si, S3, 
Vi, and V 3 couplings from — Br(n — > ev)/Br(jr — > fiv) are also shown. 



First, extend the Standard Model to include a non- 
degenerate SU(2)l triplet leptoquark (S3 or V3 in the 
notation of [55], with masses in the 0.5-1.5 TeV range. 
At the LHC these leptoquarks will be produced primar- 
ily in pairs through gluon fusion, and each leptoquark 
will decay to a lepton and a jet [72]. The peak in the 
lepton-jet invariant mass distribution will be easily de- 
tected over background. This will provide the leptoquark 
masses but yield little information about their couplings 
to fermions. The leptoquarks will also shift the neutrino- 
nucleon effective coupling g 2 L in a way that depends sensi- 
tively on both the leptoquark couplings and masses. Such 
a lcptoquark-induced shift could provide an explanation 
for the NuTeV anomaly [6TT \67\ [73]. In this scenario, 
NuSOnG would find that isospin and the strange sea can 
be constrained to the point that they do not provide an 
explanation for the NuTeV anomaly, thus the NuTeV 
anomaly is the result of new physics. The NuSOnG PW 



measurement of sin 29w will agree with NuTeV; gjj and 
the ve and Ve elastic scattering measurements will agree 
with LEP. Fig. 15 illustrates this example. NuSOnG 's 
measurement of g\ would provide a sensitive measure- 
ment of the leptoquark couplings when combined with 
the LHC mass measurements as inputs. 

A second example is the existence of a fourth genera- 
tion family. A fourth family with non-degenerate masses 
(i.e. isospin violating) is allowed within the LEP/SLD 
constraints |74j . As a model, we choose a fourth fam- 
ily with mass splitting on the order of ~ 75 GeV and 
a 300 GeV Higgs. This is consistent with LEP at la 
and perfectly consistent with Mw, describing the point 
(0.2,0.19) on the ST plot. In this scenario, LHC will mea- 
sure the Higgs mass from the highly enhanced H — > ZZ 
decay. An array of exotic decays which will be difficult to 
fully reconstruct, such as production of 6 W's and 2 b's, 
will be observed at low rates. In this scenario, isospin 
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FIG. 14: Shifts in g\ and g\ due to leptoquarks. Horizontal 
lines indicate the projected la limits of NuSOnG. 
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FIG. 15: NuSOnG expectation in the case of a Tev-scale 
triplet leptoquark. For clarity, this plot and the two follow- 
ing cases, show the expectation from only the two highest 
precision measurements from NuSOnG: g\ and v ES. 
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FIG. 16: NuSOnG expectation if the NuTeV anomaly is due 
to isospin violation and there is a heavy 4th generation with 
isospin violation. 
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FIG. 17: If LHC sees a Standard Model Higgs and no ev- 
idence of new physics, NuSOnG may reveal new physics in 
the neutrino sector. 



demonstrated. 

Lastly, while it seems unlikely, it is possible that LHC 
will observe a Standard Model Higgs and no signatures 
of new physics. If this is the case, it is still possible 
for NuSOnG to add valuable clues to new physics. This 
is because the experiment is uniquely sensitive to the 
neutrino sector. If a situation such as is illustrated on 
Fig. [17] arose, the only explanation would be new physics 
unique to neutrino interactions. 



violation explains the NuTeV anomaly, thus the NuTeV 
PW and the NuSOnG PW measurements agree with the 
i^eES measurements. These three precision neutrino re- 
sults, all with "LEP-size" errors, can be combined and 
will intersect the one-sigma edge of the LEP measure- 
ments. Fig. [16] illustrates this example. From this, the 
source, a fourth generation with isospin violation, can be 



VI. SUMMARY AND CONCLUSIONS 

NuSOnG is an experiment which can search for new 
physics from keV through TeV energy scales, as well as 
make interesting QCD measurements. This article has 
focussed mainly on the Terascale physics which can be 
accessed through this new high energy, high statistics 
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neutrino scattering experiment. The case has been made 
that this new neutrino experiment would be a valuable 
addition to the presently planned suite of experiments 
with Terascale reach. 

The NuSOnG experiment design draws on the heritage 
of the CHARM II and CCFR/NuTeV experiments. A 
high energy, flavor-pure neutrino flux is produced us- 
ing 800 GeV protons from the Tevatron. The detec- 
tor consists of four modules, each composed of a finely- 
segmented glass-target (SiG^) calorimeter followed by 
a muon spectrometer. In its five-year data acquisition 
period, this experiment will record almost one hundred 
thousand neutrino-electron elastic scatters and hundreds 
of millions of deep inelastic scattering events, exceeding 
the current world data sample by more than an order 
of magnitude. This experiment can address concerns re- 
lated to model systematics of electroweak measurements 
in neutrino-quark scattering by direct constraints using 
in-situ structure function measurements. 

NuSOnG will be unique among present and planned 
experiments for its ability to probe neutrino couplings 
to Beyond Standard Model physics. This experiment of- 
fers four distinct and complementary probes of S and T. 
Two are of high precision with the proposed run-plan, 
and the precision of the other two would be improved by 
a follow-up five-year antincutrino run. Neutrino- lepton 
non-standard interactions can be probed with an order 
of magnitude improvement in the measured effective cou- 
plings. Neutrino-quark non-standard interactions can be 
probed by an improvement in the measured neutrino- 
quark effective couplings of a factor of two or better. The 



experiment is sensitive to new physics up to energy scales 
~ 5 TeV at 95% CL. The measurements are sensitive to 
universality of the couplings and an improvement in the 
e-family of 30% and /i-family of 75% will allow for probes 
of neutrissimos. As a unique contribution, NuSOnG mea- 
sures gn/ §l, which is not accessible by other near-future 
experiments. This article described NuSOnG 's physics 
contribution under several specific models. These in- 
cluded models of Z's, extended Higgs models, leptoquark 
models and i?-parity violating SUSY models. We also 
considered how, once data are taken at LHC and Nu- 
SOnG, the underlying physics can be extracted. The 
opportunity for direct searches related to these indirect 
electroweak searches was also described. The conclusion 
of our analysis is that a new neutrino experiment, such 
as NuSOnG, would substantially enhance the presently 
planned Terascale program. 
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